Current induced vortex superlattices in nanomagnets 
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Influence of tfie spin-transfer torque on tfie vortex state magnetic nanodisk is studied numerically 
via Slonczewski-Berger mechanism. The existence of a critical current is determined for the case 
of same-directed electrical current, its spin polarization and polarity of the vortex. The critical 
current separates two regimes: (i) deformed but static vortex state and (ii) essentially dynamic 
state under which the spatio-temporal periodic structures can appear. The structure is a stable 
vortex- ant ivortex lattice. Symmetry of the lattice depends on the applied current value and for high 
currents (close to saturation) only square lattices are observed. General relations for sizes of the 
stable lattice is obtained analytically. 

PACS numbers: 75.10.Hk, 75.40.Mg, 05.45.-a, 85.75.-d 



Vortices are ubiquitous in Nature. They play a funda- 
mental role in various physical systems and one of their 
fascinating features is a possibility to create a variety of 
periodic structures. Examples are fluids [1 , supercon- 
ductors of type II where vortices and vortex lattices are 
created under the action of magnetic field [2 ; vortices 
and their arrays nucleate in rotating superfluid helium 
[3] and BEG |4j. Submicron size magnetic nanoparti- 
cles give another example where the vortices play crucial 
role. These nanoparticles due to interplay between short- 
range exchange interaction and long-range dipole interac- 
tion have a curling ground state [5 . One way to control 
the behavior of nanomagnets is to pass a spin-polarized 
current through the nanodot. As was discovered by Slon- 
czewski [6 and Berger [7 , the spin-polarized current acts 
as an effective spin-torque, which is an important feature 
for design of electrically controlled devices of spintronics 

0. 

It was shown theoretically [9UTT] and experimentally 
pT] that the dc spin-polarized current passing perpen- 
dicularly to the nanodisk can excite the circular motion 
of the vortex core. Such a circular motion can be excited 
[9 under two conditions: the current density \j\ exceeds 
some critical value jc and jap < 0, where a = ±1 is 
direction of spin polarization of the current (along the 
normal to the nanodisk £), and p = ±1 is the vortex 
polarity (the direction of the vortex core magnetization). 
When jap > 0, the vortex core remains at the center of 
the disk. 

The aim of this Letter is to show that passing the spin- 
polarized electrical current with jap > through mag- 
netic nanodisks can lead to creation of periodic vortex- 
antivortex arrays. This is a phenomenon with a thresh- 
old: it occurs when the current density exceeds some crit- 
ical value which depends on material parameters of the 
nanomagnet and its sizes. In the Letter we report about 
numerical investigation of influence of the spin-polarized 
current on vortex state nanodisk for wide range of current 
densities. Our micromagnetic simulations [12 are based 



on the Landau-Lifshitz-Slonczewski equation j6l [7t [T3]: 

rh = m X SE/Sm -\- a[m x rh] — ajem x [m x z]. (1) 

Here m = (sin cos 0, sin sin 0, cos 0) is a normalized 
magnetization, the overdot indicates derivative with re- 
spect to the rescaled time in units of {4:7T'yMs)~'^ ^ 7 
is gyromagnetic ratio, Mg is the saturation magnetiza- 
tion, and £ = E / {A-kM'^) is normalized magnetic en- 
ergy. The normalized electrical current density j = J/ Jo, 
where Jo = M^\e\h/h with e being electron charge, 
h being the sample thickness and h being Planck con- 
stant. Here and below it is assumed that the current flow 
and its spin-polarization are directed along i-axis. The 
spin-transfer torque efficiency function e has the form 
e = r]K^/ [(A^ + 1) + (A^ - l)a(m • z)] , where r] is the 
degree of spin polarization and parameter A describes 
the mismatch between spacer and ferromagnet resistance 

Studying response of the vortex state to the spin- 
transfer torque one can distinguish two critical current 
densities: ji and j2- When j < ji the stationary state of 
the system is a deformed vortex state (see below). When 
j > j2 the system goes in a saturated state when all spins 
are aligned along z axis. When ji < j < j2 a rich variety 
of dynamic states is observed: the system demonstrates 
either chaotic dynamics of vortex-antivortex 2D gas or 
regular and stable vortex-antivortex lattices. 

In the no-driving case the ground state of the nanodot 
is a vortex state with cos^^ = pfir) where an exponen- 
tially localized function /(r) describes the vortex core 
profile, and 0^ = x i 7r/2 ( (r, x) are polar coordinates 
originated in the disk center). First of all we studied how 
the ground vortex state Oy^ (j)y is modified under an in- 
fluence of the spin-polarized current. Starting from the 
current j = we increased j step by step with full relax- 
ation on each step and with size of the step Aj <C ji . A 
typical deformed vortex state is shown in the lower inset 
of the Fig. [1] Ghanges of the out-of-plane component 
of the vortex distribution are negligibly small {0 ~ Oy), 
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FIG. 1. (Color online.) Dependence of critical current Ji 
on the nanodisk thickness for different disk diameters. Insets 
show magnetization structures in nanodisk with D = 250 nm 
and /i = 15 nm for different current densities: J = 5.75 • 
10^^ A/m'^ < Ji (bottom right inset), J = 6-10^^ A/m'^ > Ji 
(upper left inset). 

while the in-plane structure is deformed appreciably and 
can be described as (/) = (/)^ + ?/^(r). The deformation 
is a non-monotonic function of r. Its amplitude is an in- 
creasing function of the current strength j. The presence 
of the deformation function ip{r) means an appearance 
of volume magnetostatic charges and increase of corre- 
sponding magnetostatic energy. As a result at j = ji, 
the vortex at the center of disk ceases to be a stable 
stationary state. The vortex escapes from the disk cen- 
ter moving in a spiral trajectory. The critical current 
ji is an increasing function of the disk thickness ( see 
Fig. [T]). This fact emphasizes the important role of the 
volume magnetostatic charges in the phenomenon. The 
mentioned spiral vortex motion is accompanied by inten- 
sive generation of the magnons and results in appearance 
of a gas of vortex- ant i vortex pairs. As a rule the mo- 
tion vortices in this state is quite irregular. However, for 
some values of the current in the range ( ji , j2 ) instead of 
chaotic gas-like dynamics stable spatially regular vortex- 
antivortex structures appear (as an example, see upper 
inset of the Fig. [T]). The following part of the Letter is 
devoted to studying this phenomenon. 

To study the observed superstructures in detail we 
used two methods: (i) the spin current of the certain 
density is sharply applied to the vortex state nanodisk; 
(ii) the nanodisk is previously saturated to the uniform 
state by a strong external magnetic field which is ap- 
plied perpendicularly to the disk plane, then the cur- 
rent of necessary density is switched on and the ex- 
ternal field is adiabatically diminished down to zero: 
B = e~^^^[BQZ + &(r',t)][l5], where Bq is saturation 
field, time parameter r is much greater than the typi- 



cal time of spin waves dissipation and additional field h 
with \h\ <^ Bq was added to avoid instable equilibrium 
states. 

Possible vort ex- ant i vortex structures which appear for 
ji < j < j2 are classified and shown in the Fig. |2j The 
obtained current diagram has the following properties: 
(i) the regular structures appear in the vicinity of the 
critical currents ji and ^2; (ii) the regions where the reg- 
ular structures exist, become wider with disk thickness 
increasing and these regions vanish for thin disks; (iii) 
regular structures have a well pronounced symmetry: a 
system with rotary- reflect ion (6*2) symmetry and systems 
with rotational Cn {n = 3, 4, 6) symmetry were observed. 
They rotate bodily around the disk center with frequency 
much lower then the corresponding gyrofrequency of the 
free vortex. The magnitude of the rotation frequency 
and its sign depend on the current; (iv) the structures 
become more compact with the current increasing; (v) 
for the currents close to j2 a crystal-like structure with a 
square lattice appears; (vi) as the current decreases line 
defects appear in such vort ex- ant i vort ex crystals, see in- 
set "f of the Fig. [2] (the defects are surrounded by dashed 
circles). As a result the crystal loses the long-range sym- 
metry but a short-range order is preserved: each vortex 
has four nearest neighbors-antivortices and vice versa. 
The defects permanently move within the structure and 
the whole dynamics is similar to the dynamics in fluids; 
(vii) further decrease of the current destroys the short- 
range order in the system and a vortex-antivortex gas 
phase appears. 

It is important to note that the crystal-like struc- 
tures appear in current driven nanodots even in the no- 
damping case {a = 0) and the interval of their existence 
only weakly depends on the damping coefficient a and 
does not depend on the method of structure creation. In 
contrast to this, the thresholds and the intervals of exis- 
tence of regular structures with a few vortex-antivortex 
pairs depend both on the damping coefficient and on the 
method of generation. These regular structures exist also 
in the no-damping case. 

The total vorticity q = Xl^i Qk^ where N is the total 
number of the particles: vortices (the winding number 
— +1), antivortices (g^^ = —1), is the same for all 
regular structures and it is equal to +1 (as for the sin- 
gle vortex). The structures with a few vortex-antivortex 
pairs also have g = 1 if one takes into account edge soli- 
tons. The edge soliton is a half of an antivortex which 
moves on the disk edge and it has the winding number 
g/c = —1/2. Each of above-mentioned regular vortex- 
antivortex structures represents a stationary state of the 
system. It means that the energy £ and the total mo- 
mentum [16 d = J {cosO — 1)(1 — d^(j))dr must be time 
independent. Considering as the most intriguing the no- 
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FIG. 2. (Color online.) The current diagrams of periodic structures in the nanodisk with diameter D — 350 nm for thicknesses 
/i = 5, 10, 15 and 20 nm. The data are obtained using numerical simulations by the method (ii), see text. For /i = 20 nm three 
additional diagrams were obtained: 1 - method (i) is used, 2 - method (ii) with zero damping (a = 0), 3 - method (i) with 
a — 0. The obtained states are classified in the following groups: "v" - static deformed vortex state, "g" - regime of the chaotic 
dynamics of the vortex- ant ivort ex gas, "f - fluid-like dynamics (see text), "s" - saturation along the current direction. The 
other states are illustrated in the insets: arrows show the in-plane magnetization distribution and the out-of-plane component 
is shown by color (gray scale). Right part of the inset "C4" and inset "f present the vortex- anti vortex lattices: rtices positions 
are shown by circles and ant ivort ices - by diamonds. All presented insets are obtained for thickness /i = 20 nm and for the 
following currents (in units 10^^ A/w?): "C2" - 5.5, "ca" - 7, "r" - 8, "C4" - 36, "c4(av)" - 5.5 for a = 0, "f - 24. 



damping case, one can obtain from equations ([T]) that 

£ = —aj j e sin^ ^dr, 

3 = (^jJ e sm^O{l-d^(l))dr (2) 

To describe the multiple vortex-antivortex structure the 
following ansatz can be used: 

N N 

cosO = ^ f{\z - Zk\), ^ = ^ (7fcarg(2; - Zfc), (3) 

k=l fc=l 

where z = re^^ and Zk = Rk^^^^ with {Rk^ ^k) being 
polar coordinates of the k-ih. particle. Substituting (|3| 
into the expressions ([2| and assuming that total area of 
the particles cores is much smaller than the disk area S 
one can obtain 

a ^ -aj^a, d-nhY, QkRl + ^h{i - ^ Qk) 

k k 

t^-TThaj^Y^QkRl^k. (4) 

k 

where ^ = 77 A^/ (A^ + 1). Thus the stationarity conditions 



£ = and a = can be satisfied if 

AT 

<7 = X] = 1, ^ QkRk = 0^ i>fe = ^ = const. (5) 

k 

Moreover, in process of the structure formation 3 = 
ao exp(— t/to), where to = is typical time of 

the formation, and systems with total vorticity q = 1 
have linear dependence between energy and momentum: 
E = Q3- For all regular structures obtained in the sim- 
ulations ^j^Qk = 1 (taking into account edge solitons), 
= const and \Yjk ^k{Rk/ ^ 10~^ with L being 
the disk radius. 

To validate our analytical considerations we carried out 
a direct stability check. To this end we first applied the 
field pulse in disk plane (red arrow), which was spatially 
localized nearby the central vortex in a circular range of 
50 nm in diameter, keeping the rest of the magnetiza- 
tion distribution artificially fixed during the field pulse 
time, see Fig. Such a field pulse shifts a vortex per- 
pendicular to the field direction. After the pulse was 
switched off the spin-current induced magnetization dy- 
namics was simulated as usual. We found that the vortex 
was turned back to its origin in a short time (about 60ps) 
and a whole structure continued its rotation, see Fig.js})). 
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FIG. 3. (Color online.) Stability of the periodic structure 
in nanodisk with D = 250 nm and /i = 20 nm. Upper row 
illustrates the stability of the central vortex, and lower row — 
the antivortex. Left column: the fixed (dark) and free (bright) 
magnetization distribution. The field pulse (red arrow) with 
the amplitude increasing up to 150 mT was applied inside a 
free area during t = 0.75 ns in order to shift a) vortex or c) 
antivortex. Right column: trajectories of the shifted vortex 
(antivortex) and its neighbors. The black arrow indicates the 
rotating direction of the structure. 



The similar picture takes place for the antivortex stabil- 
ity, cf. Figs, [s]^) and Fig. [sji). It allows to conclude that 
the superstructure is locally stable: only the strong field 
pulse can destroy the structure. 

In conclusion, we observed for the first time and stud- 
ied numerically the spatially periodic vortex-antivortex 
structures which appear in nanomagnets under influence 
of spin-transfer torque. The structure formation is only 
weekly influenced by Gilbert damping [17 , it results from 
effective dissipation induced by the spin current. We 
demonstrated analytically and proved numerically fulfill- 
ment of conditions ([5| for the periodic structures. 

Authors acknowledge computing cluster of Kyiv Na- 
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lations were performed. 
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